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Here, we present the results of theoretical analysis of the de Haas- van Alphen oscillations in quasi- 
two-dimensional normal metals. We had been studying effects of the Fermi surface (FS) shape on 
these oscillations. It was shown that the effects could be revealed and well pronounced when the 
FS curvature becomes zero at cross-sections with extremal cross-sectional areas. In this case both 
shape and amplitude of the oscillations could be significantly changed. Also, we analyze the effect of 
the FS local geometry on the angular dependencies of the oscillation amplitudes when the magnetic 
field is tilted away from the FS symmetry axis by the angle 9. We show that a peak appears at 
9 K Q whose height could be of the same order as the maximum at the Yamaji angle. This peak 
emerges when the FS includes zero curvature cross-sections of extremal areas. Such maximum was 
observed in experiments on the a — {BETS)4TIHg{SeCN)4. The obtained results could be applied 
to organic metals and other quasi-two-dimensional compounds. 

PACS numbers: 72.15.Gd,71.18.-|-y 
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I. INTRODUCTION 

Magnetic quantum oscillations are well known as a 
powerful tool repeatedly used in studies of electronic 
properties of various conventional metals [ij. Theory 
of quantum oscillatory phenomena such as de Haas-van 
Alphen oscillations in the magnetization and Shubnikov- 
de Haas oscillations in the magnetoresistivity of conven- 
tional three-dimensional metals, was derived by Lifshitz 
and Kosevich (LK) in their well-known work This 
theory was succesfully employed to extract valuable infor- 
mations concerning electron band-structure parameters 
from experimentally measured magnetic quantum oscil- 
lations. 

In the last two decades quasi-two-dimensional (Q2D) 
materials with metallic-type conductivity (e.g. organic 
metals, intercalated compounds and some other) have at- 
tracted a substantial interest, and extensive efforts were 
applied to study their electron characteristics. Magnetic 
quantum oscillations are frequently used as a tool in these 
studies 0,13. A theory of magnetic oscillations in Q2D 
materials was proposed in several works (see e.g. Refs. 
[i, S, 0, i, S [S El, El)- Significant progress is already 
made in developing the theory but there still remain some 
significant points not taken into account so far. The pur- 
pose of the present work is to contribute to the theory 
of de Haas- van Alphen oscillations in Q2D conductors by 
analyzing one of these points, namely, the effect of the 
Fermi surface (FS) curvature on the amplitude and shape 
of the oscillations. 

It is already shown that local geometrical features of 
the FS may strongly affect quantum oscillations of elastic 
constants in both conventional and Q2D metals 13l. Il4l|. 



could be easily given an explanation. In general, quan- 
tum oscillations are specified with contributions from 
vicinities of effective cross-sections of the FS. Those are 
cross-sections with minimum and maximum sectional ar- 
eas. When the FS curvature becomes zero at an effective 
cross-section, the number of electrons associated with the 
latter increases, and their response enhances. This may 
significantly strengthen the oscillations originating from 
such cross-section and change their shape and phase. 



II. THE MODEL 

The FSs of Q2D metals are known to be systems of 
weakly rippled cylinders. Accordingly, the current theory 
adopts the tight-binding approximation to describe the 
energy-momentum relation for the charge carriers. So, 
the energy spectrum could be written out as follows: 
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so one has grounds to expect similar manifestations in 
the magnetization oscillations. The effect of the FS local 
geometry on quantum oscillations in various observables 



Here, z axis is assumed to be perpendicular to the con- 
ducting layers; Px,Py,Pz are the quasimomentum p co- 
ordinates, and d is the interlayer distance. The depen- 
dence of the energies En of Px,Py is introduced in the 
Eq. (1) to take into account the anisotropies of the en- 
ergy spectrum in the conducting layers planes. Due to 
the weakness of the interlayer interactions in Q2D con- 
ductors, the coefficient Eo{px,Py) is much greater than 
the remaining coefficients in the Fourier series in the Eq. 
(1). Usually, in studies of magnetic quantum oscillations 
in Q2D conductors the general energy-momentum rela- 
tion (1) is simplified. The anisotropy of the energy spec- 
trum in the layers planes is neglected, and only the first 
two perm in the sum over "n" are kept. Then we arrive 
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FIG. 1: Schematic plots of a Fermi surface with a cosine warp- 
ing corresponding to the model (2) (left), and a Fermi surface 
of a complex profile including cross-sections with maximum 
areas where the FS curvature becomes zero (right). 



at the following simple model for a Q2D Fermi surface: 
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where p_L is the quasimomentum projection on the layer 
plane, and m± is the effective mass corresponding to the 
motion of quasiparticles in this plane. The parameter 
t in this expression (2) is the interlayer transfer inte- 
gral whose value determines how much the FS is warped. 
When t goes to zero the FS becomes perfectly cylindri- 
cal. 

Within this commonly used approximation the FS 
profile is completely established. When a strong mag- 
netic field B is applied along the normal to the layers 
(B = (0, 0, -B)) , the FS Gaussian curvature at the effec- 
tive cross-sections K^x is given by: 
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where {dP A/ dpi) ex means that the second derivative is 
calculated at the effective cross-section; A^x is the cross- 
sectional area. So, the curvature takes on values propor- 
tional to those of the interlayer transfer integral i, and 
it becomes zero only when the latter turns zero, and the 
FS passes into a perfect unwarped cylinder. 

However, there are grounds to believe that the FSs 
of some realistic Q2D conductors possess more complex 



geometries than those described by the Eq. (2). For 
instance in the experiments on magnetic quantum oscil- 
lations in the layered perovskite oxide Sr2RuOi 
it was shown that to adequately describe the FS of this 
material one must take into consideration at least four 
terms in the expansion (1). Here, we consider the Fermi 
cylinder corrugation of an arbitrary shape provided that 
the cylinder remains a surface of revolution. Separating 
out the first term in the expansion (I) we may rewrite 
the energy momentum relation in the form [l3|: 
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with e„ = —En/2t. ft follows from this expression that 
e{pzd/h) is an even periodic function of pz whose pe- 
riod equals 2ttTi/ d. Omitting all terms with n > 1 and 
assuming E — 2t, we may reduce our energy momen- 
tum relation (4) to the simple from (2). By introducing 
this expression we get opportunities to describe Q2D FSs 
of various profiles (see Fig. 1) and to analyze the influ- 
ence of their fine geometrical features on the de Haas- 
van Alphen oscillations. As shown below, these studies 
bring some nontrivial results which could not be obtained 
within the simple approximation (2). 



III. QUANTUM OSCILLATIONS IN 
MAGNETIZATION 



We start from the standard expression for the longitu- 
dinal magnetization: 



M|| (B, T, ^x) = M,{B, T,^l)^-{ — 



dB 



(6) 



Here, the magnetization depends on the temperature T 
and the chemical potential /i. The latter itself is a func- 
tion of the magnetic field and temperature, and oscillates 
in strong magnetic fields. The expression for the thermo- 
dynamic potential could be written out in a usual fashion: 

f](S, T, /i) - -T ^ In { f + exp [(/i - E)/T] } (7) 

where the summation is carried over all possible states of 
quasiparticles. When a strong magnetic field is applied 
the quasiparticles have the Landau energy spectrum, so 
the expression (4) takes the form: 
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Here, oj is the cyclotron frequency, ojq = f3B; (3 is the 
Bohr magneton, a is the spin quantum number and g 
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is the spin splitting coefficient (g -factor). Accordingly, 
we rewrite Eq. (7) as follows: 
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In further consideration we assume as usually, that the 
cyclotron quantum Hco is small compared to the chemical 
potential /i. Then we employ the Poisson summation 
formula: 

= / /W[l + 2ReE(-ir 

n=0 ^ r=l 

X exp{2Trirx)\dx. (10) 

Using this formula, the expression for the thermodynamic 
potential could be presented as a sum of a monotonous 
term fio and an oscillating correction Ail : 
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The function I{Ecr) is given by: 
I{E,) = 2iJ exp 



1 + cxp[(^, - Ai)/T] ■ 
(11) 
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where A is the magnetic length, and A{E^,pz) is the 
cross-sectional area. 

Until this point we followed LK theory in derivation 
of the expression for AQ. As a result we arrived at Eqs. 
(11), (12) which are valid for conventional 3D metals as 
well as for Q2D and perfectly 2D conductors. Diversities 
in the expressions for Afl appear in the course of calcu- 
lations of the function I{Ea). These calculations bring 
different results for different FS geometries. In deriving 
the standard LK formula it is assumed that the FS cur- 
vature is nonzero at the effective cross-sections with the 
extremal areas, and I{E„,pz) is approximated using the 
stationary phase method. For 2D metals the calculations 
of I{E„) are trivial for the FS is a cylinder and the cross- 
sectional area A does not depend of Pz- Obviously, in 
this case the FS curvature is everywhere zero. The oscil- 
lating part of the thermodynamic potential Q, for a 2D 
metal takes on the form: 



An = Nhuj 



^ [nrf 



(13) 

Here, F — cA/2TThe; N is the density of charge car- 
riers, and RT{r),Rs{r) describe the effects of temper- 
ature, and spin splitting. Also, the scattering of elec- 
trons deteriorates magnetic quantum oscillations for it 
causes energy levels broadening. The simplest way to 



account for the effects of electron scattering on the os- 
cillation amplitudes is to introduce an extra damping 
factor Rd{i") (Dingle factor) into the expression (13). 
The usual approximation for the latter has the form 
RD{r) = exp(-27rr/u;r) fl| where t is the scattering 
time of electrons. In further calculations we adopt this 
simple form for i?£)(r) for more sophisticated expres- 
sions are irrelevant to the main point of our subject. As 
a result we obtain: 
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Correspondingly, we arrive at the following result for the 
oscillating part of the longitudinal magnetization [loj : 
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(15) 



Taking into account the pz depending term in the charge 
carriers energy spectrum within the approximation (2) 
one can expand the integrand in the Eq.(12) in Bessel 
functions and easily carry out integration over pz ■ Then 
AAf|| takes the form [l3|: 



AM,, = -2NP^ y tlLRir)Jo ( ^] sin ( 2nr^ 



\ hu J V ' ^ 

(16) 

When the FS warping is negligible {t <^ huj) this ex- 
pression passes into the previous formula (15) describing 
the magnetization of a 2D metal. In the opposite limit 
(t > huj) one can use the corresponding asymptotic for 
the Bessel functions. As a result the expression for AM|| 
is transformed to the form similar to the LK result for 
the conventional 3D metals: 
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(17) 



where -Fmax — EAt^^^-^/A^ -^min — E A^j^[t^ / A, v4.j^ax; ^min 

are the maximum and minimum cross-sectional areas, re- 
spectively. As before, A is the cross-sectional area of the 
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incorrugated cylindrical FS. Before we proceed we re- 
mark again that the commonly used approximation (2) 
is not suitable to analyze the effects of the FS shape in 
Q2D metals. Within this model, the resulting formulas 
for AM are either reduced to the 2D limit [Eq. (15)] or 
they describe only cosine warped profiles of the FS with- 
out any variations [Eqs. (16), (17)]. In both cases some 
important features of the oscillations could be missed. 

We may expect the effect of the FS curvature on the 
magnetization oscillations to appeare when the FS warp- 
ing is distinct [t > Ulo). To analyze these effects we re- 
turn back to our generalized energy-momentum relation 
(4) . Then we assume that the FS curvature becomes zero 
at the effective cross-section at pz = p*. Then, as follows 
from the expression for the FS curvature (3), cPA/dpl 
equals zero at pz — p* , so we can write the following 
approximation for the cross-sectional area: 



A{pz) « Aea: ± 



(20 



dpi 



(18) 



Pz=P' 



where Z > 1. 

Performing integration in parts in the expression (11) 
we obtain 



An = 



flLO 



n^hX^ ^ (7 

r— 1 ^ 



xnjd 



R{r) I cos ( ^^A{pz) ) dpz 

(19) 

Then, using the approximation (18), and applying the 
stationary phase method to compute the integral over 
Pz in the expression (19), we obtain the following result 
for the contribution SQ from the nearly cylindrical cross- 
section to the oscillating part of fl : 
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Here, p=l + 1/2/ 
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Vo is the FS volume in a single Brillouin zone, and T{p) 
is the gamma function. Basing on this expression (20) 
we get the corresponding term in the oscillating part of 
the magnetization: 

1/21 
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(22) 



where 



1/21 



(23) 



To arrive at the complete expression for AAf|| we must 
sum up terms originating from all effective cross-sections 
of the FS. 

The FS shape near pz — p* is determined by the shape 
parameter I. When I — 1 the FS has a nonzero cur- 
vature at the considered cross-section. In this case Eq. 
(22) agrees with the well-known LK result. Assuming 
that there are two extremal cross-sections (with the min- 
imum and maximum cross-sectional areas, respectively), 
and summing up the contribution from the both we can 
easily transform our result (22) to the form (17). 

When Z 3> 1 we may roughly estimate 
\d^^e/dpf\ ^ . as {21)^^+^/(21 + 1) and (20! as 
exp[(2/ -h 2)'ln(2/ + 2)] (see [ll). So, lim,^oo 6 = 1/2, 
and Eq. (22) passes into the expression (15) describing 
magnetization oscillations in 2D conductors multiplied 
by 1/2. This extra factor appears because the ex- 
pression (22) describes the contribution from a single 
nearly cylindrical cross-section of the FS. When the 
shape parameter for both effective cross-sections goes 
to infinity, their contributions to the magnetization 
oscillations become identical, and putting them together 
we arrive at the expression (15). In general, one may 
treat I as a phenomenological parameter whose actual 
value could be found from experiments. The greater is 
the value of this parameter the closer is the FS to a 
cylinder near pz = p* ■ 

Oscillations in magnetization described by the expres- 
sion (22) vary in magnitude, shape and phase depending 
on the value of the shape parameter I which determines 
the FS local geometry near the effective cross-section. 
This is illustrated in the Fig. 2. As shown in this fig- 
ure, when there is the close proximity of the FS near 
Pz = p* to a cylinder, the oscillations are sawtoothed 
and resemble those occuring in 2D metals 0, [lH or orig- 
inating from cylindrical segments of the FSs in conven- 
tional 3D metals (2^, 21 1. When the FS curvature takes 
on nonzero value at Pz = p* (/ = 1) the oscillations are 
similar to those in conventional metals. Here, we em- 
phasize the difference between our result (22) and the 
expression (16). Using the latter one could easily obtain 
sawtoothed oscillations typical for 2D metals but only for 
small values of the transfer integral t {t ^Tiui) when the 
FS crimping is negligible. The present result (22) shows 
that the oscillation shape and phase may be determined 
not by the value of t itself but rather by the form of 
the function eipzd/Ti) specifying the FS profile. The 
sawtoothed oscillations in magnetization could occur at 
t ~ TiLj, when the FS curvature becomes zero at an effec- 
tive cross-section. To ease the interpretation of this point 
one may imagine a FS shaped as a step-like cylinder. The 



5 



0.1 



-0.1 





(Bq/B-1)x10'' 



(Bg/B-1)x10' 



FIG. 2: De Haas- van Alphen oscillations described by the Eq. 
(22) for / = 4 (left panel) and / — 1 (right panel). Calcula- 
tions are carried out for T = To = 0.5K, Bo = lOT, F/Bo = 
300; Td is the Dingle temperature, and Mo — 2N jSui/ujQ. 



curvature of such FS is everywhere zero, and oscillations 
from both kinds of the cross-sections (with minimum and 
maximum cross-sectional areas, respectively) should be 
similar to those in 2D metals. Nevertheless, the differ- 
ence in the cross-sectional areas (the FS crimping) could 
be well pronounced, and a beat effect could be mani- 
fested. Obviously, this effect is absent when t ^Tiuj and 
the FS warping is negligible. 

Also, it may happen that the FS curvature becomes 
zero at some effective cross-sections and remains nonzero 
at the rest of them. Then the contributions from zero 
curvature cross-sections (l > 1) would exceed in magni- 
tude those originating from the oridinary cross-sections 
{I = 1). This follows from the expression (22) where 
the factor {hcj/ty^^'' [Ulo < 1) is included. Depending 
on the value of the shape parameter I this factor takes 
on values between {hw/tY/"^ [I — 1) and 1 (^ ^ oo). 
So, when there is a close proximity of the FS to a pure 
cylinder at some extremal cross-sections, the contribu- 
tions from these cross-sections would be predominating, 
and they will determine the shape and amplitude of the 
magnetization oscillations in whole. 



IV. THE EFFECT OF THE CHEMICAL 
POTENTIAL OSCILLATIONS 

It is known that chemical potential oscillates in 
strong (quantizing) magnetic fields. These oscillations 
and the oscillations in the magnetization are closely re- 
lated and described by similar series. The chemical po- 
tential is determined by the equation: 



(24) 



When the quantizing magnetic field is applied, the charge 
carriers density gets an oscillating correction AiV(/i) : 



N = Noifi) + AiV(/i) 



(25) 



where A'o is the charge carriers density at B — 0. Pro- 
vided that the charge carriers density is fixed, this correc- 
tion is to be balanced by an oscillating term A/i which 
appears in the chemical potencial due to the magnetic 
field. These corrections are related to each other by the 
equation [!]]: 



(26) 



where dNo/diJ, = Dq is the charge carriers density of 
states at the FS in the absence of the magnetic field. So, 
we have: 



1 



(27) 



Assuming that the FS curvature becomes zero at the ef- 
fective cross-section at pz = p* , using the expression 
(20) for AO, and omitting smaller contributions from 
the remaining cross-sections of nonzero curvature we ob- 
tain: 



A/i = —Qhuj 



^ (-!)'• 
2-^ (7rr)P 

r— 1 ^ ^ 



Rir)sin[2nr^±^^]. (28) 



Here, Q = ANai{B / F)p /fiujDo, and Dq is the elec- 
tron density of states at the FS in the absence of the 
magnetic field. The dimensionless factor Q takes on 
small values of the order of {Tiuj / EpY^'^'' ■ This result 
(28) as well as the expression for the oscillating part 
of magnetization (22) do not contradict the correspond- 
ing results for 2D and 3D metals reported in earlier 
works 0, [ii. For a conventional 3D metal whose FS 
has a nonzero curvature at the effective cross-sections 
the oscillating correction A/i is very small compared to 
Ep [Afi/Ep - (huj/Ep)^^'^). Therefore it does not bring 
noticeable changes in de Haas- van Alphen oscillations. 
However, in the case of Q2D metals this correction could 
be more important 13, To analyze the effect of the 
quantum oscillations in chemical potential, one must take 
into account that the extremal cross-sectional areas A^x 
include corrections proportional to A/i, namely: 



Aex = AexiO) + 2TTm^Afi{B) 



(29) 



where Aex{0) is the cross-sectional area in the absence 
of the magnetic field. Correspondingly, the argument of 
the cosine function in the Eq. (20) must be written in 
the form: 



27rr- 



± — = 27rr 
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Tpex 

B 



(30) 



To simplify further analysis we keep only the first term 
in the expansion (28). Then we can employ the identity: 



exp 



zr0i?(l)sin(2^:^± J 
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(31) 
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where J„(rC;i?(l)) are the Bessel functions. 

Using this identity we can write the expression for 



AM|| as follows: 



AM|| = -2Nai(3 



F 



TB.cosAm^SliliU. ,32) 



Here, the coefficients Ag, Bg are given by: 
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7r(r + f) 
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^Jr^s{rCiR{l)): 



7r(r — f ) 
4? 



(34) 



These formulas (32)-(34) arc generalizations of the re- 
sults obtained for 2D metals [11, 23, 24]. 

As follows from Eqs. (32)-(34) the oscillating correc- 
tion to the chemical potential may bring some changes 
in the amplitude and shape of the de Haas-van Alphen 
oscillations in Q2D metals. Keeping in mind that the 
damping factor R{r) takes on values less than unity, 
and reduces while r increases, we can write explicit ex- 
pressions for a few first harmonics in the form: 

5Mi ^ 2NaiR{l)Ti-Pii\n(^Ti^±^- (35) 
SM2 = -2iVa,i?(2)(27r)"''<^ sin U7r-|- ± ^ j 



R{2) ^ \ B 21 



SM3 = 2A^a,i?(3)(37r)"''| sin (^67r^ ± 

i?(2)i?(l) fSW . / , TT 



(36) 



(37) 



So we see that the chemical potential oscillations do not 
affect the fundamental harmonic but they contribute to 
higher harmonics bringing changes to their amplitude 
and phase. Similar conclusions were recently made ana- 
lyzing the effect of chemi cal p otential oscillations on mag- 
netization in 2D metals [ll|. The effect depends on the 
FS local geometry. When the shape parameter I takes 
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FIG. 3: The effect of the chemical potential oscillations on 
the magnetization oscillations. The curves are plotted taking 
into account the chemical potential oscillations according Eq. 
(32) (left), and neglecting the latter (right). For both curves 
i = 4 (p = 1.125) and the remaining parameters coincide 
with those used in plotting Fig. 2. 



on greater values, the effect becomes stronger. As shown 
in the Fig. 3, for a pronounced proximity of the FS to 
a cylinder near the effective cross-section {I = A), a no- 
ticeable difference in the magnetization oscillations shape 
and magnitude arises due to the effect of the chemical po- 
tential oscillations. 

The correction from A/i brings changes in the posi- 
tion of spin-splitting zeros in the harmonics of the de 
Haas- van Alphen oscillations. These changes were stud- 
ied before within the simple approximation (2) for the 
energy spectrum [26] . It was shown that the spin-zero 
positions for the second and third harmonics are not 
completely determined with zeros in corresponding spin- 
splitting factors Rs{2), i?s(3). They also depend on the 
values of R{1),R{3) (which are dependent of tempera- 
ture and scattering), and on the magnetic field. Here, 
we show that the spin-zero positions also depend on the 
FS geometry at the extremal cross-sections. This follows 
from the expressions (36), (37). 



V. ANGULAR DEPENDENCE OF THE 
OSCILLATIONS AMPLITUDES 

The effect of the FS curvature on the quantum oscil- 
lations in the magnetization is expected to be very sen- 
sitive to the geometry of the experiments. The reason 
is that the effective FS cross-sections (with the mini- 
mum/maximum cross-sectional areas) run along the lines 
of zero curvature (if any) only at certain directions of the 
magnetic field. When the magnetic field is tilted away 
from such direction by the angle 9, the extremal cros- 
section slips from the nearly cylindrical strip on the FS 
containing a zero curvature line. This results in the de- 
crease in the oscillations amplitude. The phase of the 
oscillations also changes. These angular dependencies of 
the oscillations amplitudes and phases radically differ in 
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origin from the effect first described by Yamaji 27| . 

The Yamaji effect occurs due to the coincidence of the 
FS extremal areas Amax and A^i^ at certain angles of 
inclination of the magnetic field with respect to the FS 
symmetry axis. At such angles all the cross-sections on 
the FS have the same area, so the amplitude of the de 
Haas-van Alphen oscillations increases. The Yamaji ef- 
fect originates from the periodicity of the Pz dependent 
contribution to the charge carriers energy spectrum, and 
it is inrelated to the presence/absence of zero curvature 
lines on the relevant FS. Also, there is a crucial difference 
in the manifestations of the two effects. The angular de- 
pendence originating from the effects of the FS curvature 
reveals itself at the very small values of 9, whereas the 
first maximum due to the Yamaji effect usually appears 
a.t 9 10° or even greater. To further clarify the differ- 
ence between the two effects we analyze the angular de- 
pendence of de Haas-van Alphen oscillation amplitudes 
assuming that the FS curvature becomes zero at an ex- 
tremal cross-section when the magnetic field is directed 
along the FS axis of symmetry. 

We suppose that the magnetic field is inclined from 
the FS symmetry axis by the angle 9 within the xz 
plane, and we employ the coordinate system whose z' 
axis is directed along the magnetic field. We use the 
energy momentum relation given by Eqs. (4), (5), and we 
rewrite them in terms of new coordinates p'^,py^ p'^ {p'^ = 
Px cos 9 + Pz sin 9\ p'^ — pz cos 9 ~ px sin 9) . Changing 
variables in the Eqs. (4), (5) and keeping in mind that 
sm9 takes on very small values at small angles 0, we 
may present the energy momentum relation for a weakly 
corrugated Fermi cylinder {t fi) in the form: 



pI = Px cos^ 9 +pI- 4:m±t 



nd, , 
e„, cos I — [Pz cos ( 



E 

n=l 



p'xsm9]] (38) 



where po is the radius of the inwarped Fermi cylinder 
(po — \J Ajii). Introducing polar coordinates in the cut- 
ting plane (p^ — p cos ip; py — psunp) we may calculate 
the FS cross-sectional area provided that the magnetic 
field is tilted away from the FS symmetry axis: 

A{p'z, 9) = r I'pdp - Ao{9) + AA{p',,9). (39) 



Here, 



Ao{9) = 2pl 



ir/2 



dip 



cos^ 9 cos^ (p + sin^ tp cos ( 



(40) 



and 



AA(p^,cos0) = 2mj^i^e„y cos |^ ^ cos 9 



npod 



tan 9 cos ip ) dp 



(41) 



Then we can present the oscillating part of the longitu- 
dinal magnetization in the form: 



i-iy 



X sm 



r=l 

2TirF{9) 
B 



^r{9) Yr{9). (42) 



Here, Yr{9) = ^C^{9) + S^{9), 

$ = arctan[5^(e')/a(^)], F{9) = F/cos9, 



Cr{9) = 



2TTh COS 9 



(7Th/d}cose /y2 

COS —j^l\A{p'^, 

(-Trn/d) cose V n 

{'KTij d) cos 



dp'z, 
(43) 



SriO) = :^-^n I smr-^AA{p'^,9)\dp',, 

(44) 

Expanding the integrand in (41) in Bessel functions we 
can easily carry out the integration over tp. Then we get: 



A A{p'z, cost 



s—^ I lip, 

4:nm±t cos y — ^ cos ( 

n=l 

' np^d 



Xjr 



■ tan( 



(45) 



The first term in this expansion coincides with the re- 
sult given in the Yamaji work. The latter was obtained 
assuming the simple cosine warping of the FS described 
with the Eq. (2). 

To analyze the effect of the FS curvature we assume 
that the curvature becomes zero at pz = T:h/d. Requir- 
ing that {d^A/dp 



(see Eq. (3)) and keep- 
ing only two first terms in the expansion (5) we obtain: 



' p^—'Kh / d 



Pzd 



cos 



Pzd 



1 f2pzd 

— cos 

4 \ h 



(46) 



In this case the cross-sectional area A{pz) near pz = 
Trh/d is approximated as: 



A{pz) = A{po) + 



TTm^t / d 

n 



Pz 



7rh\ 

~d) 



(47) 



Comparing this expression with the Eq. (18) we conclude 
that the shape parameter 1 = 2. 

Correspondingly, we must put ei = 1, 62 = 1/4, e„ = 
(n > 2) into the general expression (46) for Ayl(p^, 9). 
As a result we get: 



A^(p'^,cos( 



1 



I ^^^^ a\ T f Po'^ i- , 
cos I -j^ COS w I Jo I tan ( 



2pU .\ , f2pod 

cos I COS 9 \ Jn \ tan 9 

4 V ^ J \ h 



(48) 
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FIG. 4; Angular dependencies of the magnetization oscil- 
lations amplitudes. The curves are plotted assuming that 
t/hu} = 0.5, pod/h = 47r. The shape parameter I takes on 
the values: I — 1 (solid line); / — 2 (dashed line); I — 3 
(dotted line). The plotted curves are described by the Eqs. 
(42)-(44) (r = l). 



To describe FSs possessing closer proximity to a perfect 
cylinder near the certain extremal cross-section we must 
keep more terms in the expansion (5). For instance, as- 
suming ei = 1, 62 = 2/5, £3 = 1/15, and e„ = (n > 
3) we ensure that both (PA/dpl and d^A/dp^ become 
zero at pz = Trh/d, which corresponds to 1 — 3. Simi- 
larly, at ei = 1, £2 = 1/2, £3 = 1/7, £4 = 1/56, and 
£„ = (n > 4) we obtain I — 4, and so forth. Sub- 
stituting these numbers into the general expression (45) 
for AA{p'z,cos6), and using the results to compute the 
functions Cr{6) and Sr{0) given by the Eqs. (43), (44), 
we may finally calculate the factors Yr{0). The latter de- 
scribe the desired angular dependence of the oscillations 
amplitudes. 

Here, we carried out the calculations accepting t/hcu — 
0.5 and pad/fi = Att and keeping only the first term in 
the sum over r in the Eq. (42). The resulting curves are 
presented in the Figure 4. The solid line in this figure is 
associated with the energy spectrum of the form (2). The 
corresponding FS has a cosine warping, and possesses a 
nonzero curvature at the expremal cross-sections. The 
high peak at 6* = 0.185 (10.6°) corresponds to the first 
Yamaji maximum. The position of this peak is in agree- 
ment with the equation (pod/h) tanO = 37r/4 (see [27|). 
Two preceding zeros originate from the beats. The re- 
maining curves represent FSs whose curvature becomes 
zero at their minimum cross-sections at = 0. We see 
that the closer the FS shape to a that of a perfect cylin- 
der in the vicinities of these cross-sections (the greater is 
the value of I ) the greater is the amplitude of the oscil- 
lations near 9 = and the smaller is the Yamaji maxi- 
mum. At ^ = 3 the Yamaji maximum is approximately 
2 times higher than the maximum at 6 = 0, whereas 



at ^ = 1 the ratio of the heights takes on the value 
close to 4. We may expect that at very close proximity 
of the FS to a cylinder near the extremal cross-section 
(I ^ 10), the amplitude maximum at 9 = will exceed 
the Yamaji peak. Also, we see that the amplitudes of 
the oscillations associated with the FSs containing zero 
curvature extremal cross-sections do not become zero at 
small angles 6. This is due to the fact that the consid- 
ered FSs are warped cylinders whose curvature is zero 
at the cross-sections with the minimum areas (such as 
Pz = nh/d) but remains nonzero at the cross-sections 
with the maximum areas (such as Pz — 0). Difference in 
the amplitudes of the oscillations originating from the FS 
local geometry near its extremal cross-sections prevents 
beats from being well manifested. 

The angular dependence of the magnetization oscil- 
lations amplitude resembling that presented in the Fig. 
4 was reported to be observed in experiments on the 
Q2D organic metal a - {BETS)2TIHg{SeCN)i 
In these experiments a high peak in the amplitude was 
observed when the magnetic field was directed along the 
axis of the corrugated cylinder which is the part of the 
FS. When the field was tilted away from this axis by 
the angle 9 the amplitude was rapidly reducing, and 
it reached approximately one half of the initial value 
at ~ 5°. Further increase in the angle 9 brought 
small variations in the amplitude until another peak was 
reached at 6* ^ 18°. Identifying this second peak with 
the first Yamaji maximum we may conjecture that the 
higher peak at 6* = arises due to the presence of the 
FS extremal cross-sections of zero curvature. The rela- 
tion between the heights of the peaks reported in 
gives grounds to expect that the nearly cylindrical seg- 
ments of the a-{BETS)2TIHg{SeCN)i Fermi surface 
(where the FS curvature becomes zero) are very close to 
perfect cylinders, so that the shape parameter I takes 
on values significantly greater than unity. 

VI. CONCLUSION 

In summary, the present work aims to study the effect 
of the FS shape on the de Haas-van Alphen oscillations 
in Q2D conductors. Such analysis is important for there 
exists a great deal of interest in studies of band-structure 
parameters and other electronic properties of these mate- 
rials. Usually, the simple model for the electron spectrum 
(2) is employed to extract the relevant information from 
the experiments. This approximation has its limitations, 
so some problems arise in interpreting the experimen- 
tal data (see e.g. Refs. [l^, [la, 13 )■ An important 
limitation of the current theory is that the latter misses 
the effects of the FS geometry assuming a simple cosine 
warping of the FS. Here, we lift this restriction on the FS 
shape. We show that the FS profile may significantly af- 
fect the quantum oscillations in magnetization if the FS 
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curvature becomes zero at a cross-sectional area. Also, 
we show that main characteristics of the oscillations are 
determined by two different factors, namely by the FS 
curvature at the effective cross-sections and the trans- 
fer integral, whereas existing theory takes into account 
only the latter. These two factors work simultaneously, 
and their effects could be separated. As discussed above 
the features typical for 2D conductors could be revealed 
when the FS is rippled provided that its curvature turns 
zero at the cross-sections with extremal areas. 

The presence of zero curvature effective cross-sections 
noticeably affects the angular dependencies of the oscil- 
lation amplitudes. A maximum originating due to the 
FS local geometry could emerge. The height of this 
peak could be comparable with the height of the well 
known Yamaji maximum and even exceed the latter. 
This agrees with the experimental results observed on 
the Q2D a - {BETS)2TIHg{SeCH)i organic metal 
[2^ . The proposed approach could be useful in analyzing 
experiments on magnetization oscillations in Q2D con- 
ductors, especially those where sawtooth features in the 
oscillations are well pronounced, ft could help to extract 
important extra informations concerning fine geometrical 
features of the FSs of such materials. 
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